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Abstract 

We present analyses to determine the fundamental parameters of the Galaxy based on VLBI astrometry of 52 
Galactic maser sources obtained with VERA, VLBA and EVN. We model the Galaxy's structure with a set of 
parameters including the Galaxy center distance Rq, the angular rotation velocity at the LSR £Iq, mean peculiar 
motion of the sources with respect to Galactic rotation (U stc , V SIC , W SIC ), rotation-curve shape index, and the V 
component of the Solar peculiar motions Vq. Based on a Markov chain Monte Carlo method, we find that the 
Galaxy center distance is constrained at a 5% level to be Rq = 8.05 ± 0.45 kpc, where the error bar includes both 
statistical and systematic errors. We also find that the two components of the source peculiar motion U src and W src 
are fairly small compared to the Galactic rotation velocity, being U STC = 1.0 ± 1.5 km s _1 and W src = — 1.4 ± 1.2 
km s _1 . Also, the rotation curve shape is found to be basically flat between Galacto-centric radii of 4 and 13 kpc. 
On the other hand, we find a linear relation between V SIC and Vq as V SIC = Vq — 19 (±2) km s _1 , suggesting that the 
value of V SIC is fully dependent on the adopted value of Vq. Regarding the rotation speed in the vicinity of the Sun, 
we also find a strong correlation between ilo and Vq. We find that the angular velocity of the Sun, f2 Q , which is 
defined as f2 Q = fio + Vq/Rq, can be well constrained with the best estimate of i7 Q = 31.09 ± 0.78 km s _1 kpc -1 . 
This corresponds to 0o = 238 ± 14 km s _1 if one adopts the above value of Rq and recent determination of ~12 
km s _1 . 

Key words: astrometry — VLBI — The Galaxy — Galactic parameters 



1. Introduction 

As a first-order approximation, the Milky Way Galaxy can 
be regarded as an axi-symmetric system with a circular rota- 
tion, and its structure and motions can be described combina- 
tions of fundamental parameters. For instance, historically the 
Milky Way's rotation was discovered based on the determina- 
tions of so-called Oort constants A and B, which are combi- 
nations of the distance to the Galaxy center Rq, and the rota- 
tion velocity of the LSR (Local Standard of Rest) 6o, and their 
derivatives. Today the Galactic constants Rq and Qq are most 
frequently-used parameters which provide fundamental scales 
of the Galaxy's size and rotation speed. The IAU recommenda- 
tion values are 8.5 kpc and 220 km s" 1 , (Kerr & Lynden-Bell 
1986), but accurate values of these parameters are still at issue: 
for instance, most of recent studies of Rq favor values some- 
what smaller than the IAU value (Reid 1993: Ghez et al. 2008; 



Gilssen et al. 2009; McMillan & Binney 2010), being be- 
tween 8.0 and 8.5 kpc. The rotation curve is another important 
quantity to describe the rotation and mass distributions in the 
Galaxy, and often parameterized with a few parameters (e.g., 
flat rotation curve or power-law rotation curve). However, the 
exact shape of the rotation curve also remains yet to be studied. 

Astrometry of Galactic sources is believed to be most pow- 
erful to determine the fundamental parameters of the Galaxy 
because it provides full six-dimensional phase-space informa- 
tion of Galactic objects. Most notably, until recently most 
of parallax measurements have been limited to local sources, 
within ~1 kpc from the Sun (for instance, such as using 
HIPPARCOS). However, recent development of high-accurate 
astrometry with VLBI is now providing tens of sources for 
which precise distances and proper motions are known. Thanks 
to this, new determination of fundamental parameters of the 
Milky Way Galaxy becomes feasible. For instance, Reid et 
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al. (2009b) used 18 maser sources to explore the structure of 
the Milky Way and obtained the Galactic parameters such as 
i?o = 8.4±0.6kpc and G = 254± 16km s _1 . In addition, they 
suggested that rotation of star-forming regions could lag behind 
the true Galactic rotation by ~15 km s _1 . There are also sev- 
eral studies of Galactic structures based on reanalysis of Reid et 
al.' data (e.g., Bobylev & Bajkova 2010; McMillan & Binney 
2010), which emphasize importance of Galactic scale astrom- 
etry. Since VLBI astrometry has been now extensively carried 
out with VERA and VLBA, the number of sources for which 
astrometric information is available is significantly increased 
(e.g., there are ~10 new sources reported in the VERA spe- 
cial issue in PASJ 63, 201 1). Therefore, it is worth conducting 
a new study of fundamental parameters of the Galaxy based 
on these new results. To carry out this, in this paper we com- 
pile the most-updated list of VLBI astrometry results, and pro- 
vide new analyses of fundamental parameters of the Milky Way 
Galaxy. 

2. Sample 

We collected results of VLBI astrometry from the litera- 
tures including most recent results, and compiled a list of 
parallaxes and proper motions in table 1. The list con- 
sists of star-forming regions with maser sources or contin- 
uum sources (such as T-Tau stars) that have been observed 
with VERA (VLBI Exploration of Radio Astrometry )/VLBA 
(Very Long Baseline Array )/EVN (European VLBI Network). 
The star-forming regions included here are those with H2O 
and/or CH3OH maser emissions, or gyro-synchrotron emis- 
sions which originate from magnetic activity of pre-main se- 
quence stars. We include in our sample almost all published 
works available to date: the only exceptions are: 1) Sgr B2 
(Reid et al. 2009c), which is too close to the Galaxy center (ex- 
pected to be at ~ 100 pc) and is known to have significant pecu- 
liar motion, making this sources unsuitable for the study of the 
Galactic rotation, 2) G34.43+0.24 (Kurayama et al. 2011) for 
which the proper motion in declination is not measured due to 
poor UV coverage for this nearly-0 declination source, and 3) 
IRAS 05137+3919 (Honma et al. 201 1), which is likely to be 
located in the far outer Galaxy (R > 15 kpc), where our simple 
parameterization of Galactic rotation (e.g., power-low rotation 
curve) may not be appropriate. 

We only selected star-forming regions, and excluded late- 
type stars because properties of these populations are quite dif- 
ferent from those of star-forming regions: for instance, late- 
type stars are known to have a larger velocity dispersion than 
that of star-forming regions, and their rotation lags behind 
the Galactic rotation, which is known as "asymmetric drift". 
However, we note that we included two super-giant sources 
VY CMa (Choi et al. 2008) and S Per (Asaki et al. 2010) since 
they are massive stars that evolve quickly after their birth, and 
thus can be considered as "young" populations. 

Table 1 summarizes the 52 sources that are selected in the 
manner described above. In order to show the distribution of 
the sample sources in the Galaxy, in figure 1 we show maser 
positions in the CO l-v diagram obtained by Dame et al.(2001). 
As can be seen in figure 1, the sources are mostly located in the 
northern part of the Galactic plane. This is because the ob- 



serving arrays (VERA, VLBA and EVN) are all located in the 
northern hemisphere. However, except for the concentration to 
the northern hemisphere, the sources are scattered well from 
the Galaxy center to the region beyond the anti-center, up to 
I ~ 240°. To fill up the gap in the southern hemisphere, as- 
trometry of maser sources with LBA (Australian VLBI array) 
and forthcoming SKA (Square Kilometer Array) will be effec- 
tive, and this will be one of the most important future tasks in 
this kind of study. 

As we will see in later sections, the whole set of 52 sources 
may include "outliers," which have large deviations from ex- 
pected motions. Such a source could introduce systematic error 
into our results, and thus its effect should be investigated care- 
fully. To do this, in later section (section 5), we will present 
results of various sample sets, in which some sources in list 1 
are discarded depending on its degree of deviations. The source 
elimination will be made based on source's contribution to the 
likelihood, Pj(Oj|M) (see equation [3]): we create a new set 
of sample with N — 1 sources by removing the source with 
smallest Pj(Oj|M) in the iV-source sample. In this way, we 
will try various sets of samples with source numbers ranging 
from 52 to 44 (for details see section 5). In table 1 we mark the 
sources with their ranks in the likelihood. 

3. Galaxy model 

In order to derive fundamental Galaxy parameters, we use a 
simple Galaxy model in which motions of sources are basically 
in circular rotation with small systematic/random motions. As 
is widely used, the Local Standard of Rest (LSR) is set to be in 
circular rotation around the Galaxy center, which is described 
with two Galactic constants Rq and Go. Obviously Rq is the 
distance between the Galaxy center and the LSR, and Oo is the 
Galactic rotation velocity at the LSR. The ratio of 80 to Rq 
gives the Galactic angular velocity at the LSR Qq (=Go/i?o)- 
In the following analyses, we use a set of Rq and f2o> rather 
than a set of Rq and Go, as the Galactic constants to be solved 
because it is known that Rq and Go are tightly correlated (e.g., 
Reid et al.2009b; McMillan and Binney 2010). 

Regarding the Galactic rotation curve, we adopt two differ- 
ent rotation curves, a power law model and 2nd-order polyno- 
mial model. In the power law model (which is mainly studied 
here), the Galactic rotation velocity at any radius R is assumed 
to obey the following rotation low, 

e(jR) = e °(}f )"■ (1) 

Here a is the rotation curve index, and this parameter is to be 
solved in the following analysis to determine the shape of the 
rotation curve. Obviously, a flat rotation curve with G(i?) = Go 
corresponds to a = 0. Note that R is the radius in the cylindrical 
coordinate system defined with (R, /3, Z), with its origin located 
at the Galaxy center (the notations here follow those in Reid et 
al. 2009b). On the other hand, in order to assess the effect of 
rotation curve models on Galactic parameter determinations, 
we also use the polynomial model described as, 

e(R) = Go + aQ(R-RQ) + bo(R-R ) 2 . (2) 

Here a = dQ/dR and b = (d 2 Q/dR 2 )/2, respectively, and 
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Fig. 1. Location of 52 maser sources for which accurate astrometric data are available (table 1), superposed on the longitude-velocity diagram of CO 
(Dame et al. 2001). 



both coefficients are defined at R = Rq. This rotation curve 
contains two parameters while the power law model has only 
one. Thus, the polynomial rotation curve has a potential to 
give a better fit to the data, and we can test if inclusion of an 
additional parameter would alter the results based on the power 
law model. 

We also include possible systematic motions of star-forming 
regions as model parameters to be solved, since Reid et 
al. (2009b) reported that high mass star-forming regions would 
lag behind Galactic rotation by ~ 15 km s _1 . In the present 
paper, this effect is re-considered in the same manner as Reid 
et al. (2009b) by including mean systematic motions of star- 
forming regions (U SIC , V STC , W SIC ). Note that the directions of 
Usrc, Vstc, Wsxc are defined at each source position, with U SIC 
toward the Galaxy center, V src toward the direction of Galactic 
rotation, and W 8TC toward the north Galactic pole, respectively. 

In addition to the above parameters that describe the mo- 
tions of maser sources in the Galaxy, we add one additional 
parameter to be solved for in our analysis: the V compo- 
nent of the Solar motion with respect to the LSR (V©). As 
is widely known, the Sun is not at rest in the LSR frame but 
has its own peculiar motion with respect to the LSR. Based 
on the HIPPARCOS observations of local stellar motions and 
distances, Dehnen and Binney (1998) obtained the Solar mo- 
tion with respect to the LSR as, (U Q , V Q , W Q ) = (10.00±0.36, 
5.25±0.62, 7.17±0.38) km s~\ reducing the V component of 
the Solar motion by ~ 10 km s _1 from the classical value of 
^ =15.4 km s" 1 (Kerr & Lynden-Bell 1986). However, re- 
cent reanalyses of the data reported that the V component of 
the Solar motion by Dehnen & Binney (1998) was an under- 
estimate and that an upward modification of V© is required 
by - 7 km s" 1 (Schonrich et al. 2010), yielding V Q ~ 12 
km s _1 . Interestingly, this upward modification reduces the 
rotation velocity lag of the star-forming regions (McMillan & 
Binney 2010). Therefore, the lag parameter V^ rc and V© are 
correlated to each other to some extent, and we include V© as 
possible parameter (fixed or solved depending on models con- 
sidered) to see how our results are affected. Note that the other 
two components of the Solar motion (£/©, Wq) are assumed 
to be well-determined, since previous studies were consistent 



between each other (e.g., Kerr & Lynden-Bell 1986; Dehnen & 
Binney 1998; Aumer & Binney 2009; Schonrich et al. 2010). 
Throughout the present paper, we simply take U sun and W sun 
of Dehnen & Binney (1998). Table 2 summarizes the parame- 
ters to be solved for in the present analyses. 

4. Analysis 

To obtain the best estimates of the Galactic parameters 
described above, we utilize a Markov chain Monte Carlo 
(MCMC) method with the Metropolis algorithm. To do this, 
first we define the likelihood function as 

L(M) = JJP<(0<|M). (3) 

i 

Here Pj(Oj|M) is the conditional probability of finding ob- 
servables O; given a set of model parameters M. repre- 
sents the observables for the i-th source, a detailed description 
of which is provided later. The model M represent the model 
parameters described in the previous section: for the power law 
rotation curve, 

M=(R ,n a ,a, U SIC , V SIC (or V Q ) , W src ) , (4) 

or for the polynomial rotation curve case, 

M = {Ro , O , a , b , U src , V SIC (or V©), W SIC ) . (5) 

Note that as we will see later there is a strong correlation be- 
tween V src and Vq, and thus these two parameters are not 
solved for at the same time, but one of them are fixed while 
the other is being solved for. 

In a practical MCMC method, we calculate the likelihood 
L(M) based on the equation (3) for a current set of model pa- 
rameters M. Then, the next set of M' is calculated by adding 
a small random fluctuations to the previous M, and the likeli- 
hood L'(M') is calculated. If i'(M') > L(M), the next pa- 
rameter M' is accepted. If not, we draw a random variable 
u, which has a uniform probability between and 1, and we 
accept the next parameter M' in case of L'(M')/L(M) > u. 
In other case (i.e., Z/(M')/L(M) < u), the next parameter set 
is rejected and the parameter set remains at the previous one 
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parameter notes 



Ro the distance from the Sun to the Galaxy center 

fio Angular rotation velocity at the LSR (flo = O0/-R0) 

a rotation curve index defined as <d(R) = G>o(i?/i?o) Q 

a , b rotation curve coefficients defined as 0(i?) = Oo + a (R — R ) + bo(R — Rq) 2 

U sl - C mean peculiar motion of SFRs toward the Galaxy Center 

V src mean peculiar motion of SFRs toward the direction of Galactic rotation 

JV src mean peculiar motion of SFRs toward the north Galactic pole 

Vq V component of the Solar motion (toward the Galactic rotation) 



M. Above procedures are iterated for a large number of trials 
(typically for ~ 10 5 ). The best estimates of M and their error 
bars are obtained by taking means and standard deviations of 
the samples of M in the MCMC trials, after removing early 
trials in the initial "burn-in" phase, which is set to be the first 
10% of the trials. Note that these procedures correspond to a 
MCMC method with a uniform a priori probability distribution 
of model parameters, namely P(M) = const. The reason for 
adopting a constant prior probability is that in the present paper 
we want to constrain Galactic parameters based only on VLBI 
astrometry data of star-forming regions rather than partially re- 
lying on other results in previous studies. However, since the 
parameters M are well constrained to relatively narrow ranges 
in the final results of MCMC (see next section) the choice of 
a different a priori probability (such as uniform probability in 
the logarithm) would not alter the main results presented here. 
For the observables, O; (which is for the i-th source) is given 

as 

; = (Vlsr, M> Hb)i- (6) 

Here Vlsr is the radial velocity with respect to the LSR (based 
on the conventional Solar motion), \i\ and are the two com- 
ponents of proper motions along the Galactic coordinate I and 
b. We adopted these quantities as observables because they 
are commonly used, but here we note that they are not in the 
exactly same frame: while Vlsr is measured in the conven- 
tional LSR frame, pi and pb are observed in the heliocentric 
frame. When MCMC analyses are done in later section, of 
course ui, pb and Vlsr are converted to the appropriate LSR 
frame by using the adopted Solar motion. We also note that 
VLBI astrometry provides the full 6-dimensional coordinates 
in the phase space, but here the source coordinates I and b are 
not considered as "observables". 

In order to calculate the conditional probability -P;(0;|M) 
in equation (3), here we simply adopt Gaussian for the proba- 
bility distribution of finding an observable given the model M 
and the source parallax 7T;, namely, 

rtl l»j \ 1 ( (^obs ^modcl) \ 

P( Xi \M,n) = 7 =e W { - 2 j, (7) 

where cc b s is the observable (corresponding to Vlsr, 
and pb of the z-th source). The model value x mo d c i = 
^modci(M, TTi) is calculated based on the Galactic model M 
as well as the observed source parallax m. The conditional 
probablity Pi(Oi|M) in equation (3) can be obtained by mul- 
tiplying the Gaussian distribution in equation (7) for the three 



parameters (Vlsr, W> and n b) and also by integrating over the 
possible parallax range of tt. The integration over the possible 
range of tt is done by adopting Gaussian distribution of tt with a 
standard deviation given by the observational uncertainty listed 
in table 1. Then, the likelihood in equation (3) is evaluated by 
multiplying the conditional probability Pi(Oi|M) for all the 
sources in the sample. 

The standard deviation a in equation (7) is determined by 
considering both the errors in the observations and the uncer- 
tainties in model estimates. In addition to the observational 
errors listed in table 1, we include the effect of source random 
motion, which is assumed to be isotropic with a 1 -dimensional 
velocity dispersion of a v . This velocity dispersion accounts 
for the turbulent motion in the molecular clouds and also sys- 
tematic difference of maser source motions from proto-stars, 
and throughout the paper we set a v = 7 km s _1 unless oth- 
erwise noted. As for radial velocities, observational error of 
the LSR velocity <7LSR,obs and the source velocity dispersion 
a v is summed in quadrature to calculate a in equation (7) as 
°lsr = °lsr obs + °v- F° r P ro P er motions, again we add the 
observational error er b s and the random motions in quadrature 

to obtain 07, l<h = a 2 hs ~ + (a v /D) 2 . 
5. Results 

In order to obtain the best estimates of the Galactic parame- 
ters, we conducted several MCMC analyses for different sam- 
ple data sets as well as different set of model parameters. First 
in section 5.1, we present the results with different samples 
(while fixing model parameters), and then in section 5.2 we 
treat cases with different model parameters. 

5.7. Results with varying number of sample sources 

Table 3 summarizes the MCMC results for different sam- 
ples, with the source number varying from 52 to 44. Here in 
table 3, the model parameters are Galactic constants Rq and 
r^o, the power-law index of the rotation curve, a, the mean 
SFRs motions (with respect to the Galactic rotation) U src , V src 
and W STC . Note that the Solar motion V© is fixed to 5.25 km _1 , 
adopting the results by Dehnen & Binney (1998). The effect of 
varying V© and the rotation curve shape will be examined later. 
Error bars in table 3 show statistical error (1-er) obtained with 
the MCMC analyses. In table 4, we also show the correlation 
coefficients between Galactic parameters for the sample of 49 
sources as an example (result 4 in table 3). As seen in table 4, 
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Table 3. MCMC results of the Galactic parameters for various sample sets 



result 


if or 


Kq 


S2 


a 


J J 

^src 


T / 

^src 


W src 


T T 

'© 


X /a. o.i. 


ID 


source 


(kpc) 


(km/s/kpc) 




(km s ) 


(km s ) 


(km s ) 


(km s ) 




1 


52 


o ^ 'i i n ill 

8.27±0.42 


29.98±0.64 


+0.05±0.02 


l.lzbl.3 


— 14.8±1.2 


— 0.6±1.0 


5.25* 


494/150 


9 




c 70+0 4.1 


Z7.07inU.UU 


_i_n 0S+0 0? 

T U.UJXU.Ui 


W. 0_J_ 1 . i 


It.UX 1 .z, 






4.10/14.7 

T - 1 \J/ It/ 


3 


50 


8.40±0.42 


30.37±0.76 


+0.04±0.02 


0.4±1.4 


-14.3±1.2 


-2.0±1.1 


5.25* 


388/144 


4 


49 


7.84±0.40 


30.61±0.75 


+0.00±0.02 


0.6±1.4 


-13.1±1.2 


-1.9±1.0 


5.25* 


353/141 


5 


48 


7.90±0.40 


30.33±0.71 


+0.00±0.02 


1.4±1.4 


-13.6±1.2 


-1.8±1.1 


5.25* 


333/138 


6 


47 


7.91±0.38 


30.73±0.76 


+0.00±0.02 


1.4±1.4 


-13.6±1.3 


-1.5±1.1 


5.25* 


327/135 


7 


46 


7.94±0.42 


30.59±0.76 


+0.01±0.02 


1.2±1.4 


-13.2±1.3 


-1.3±1.1 


5.25* 


312/132 


8 


45 


7.93±0.42 


30.70±0.77 


+0.02±0.02 


0.8±1.5 


-12.7±1.3 


-l.Oil.l 


5.25* 


303/129 


9 


44 


8.04±0.42 


30.71±0.76 


+0.02±0.02 


0.8±1.5 


-12.5±1.3 


-0.8±1.1 


5.25* 


288/126 



: fixed parameter 



For result 1, all 52 sources are used. For other samples, following sources are removed one43y-one from result 2 to 9: IRAS 
20126+4104, G48.61+0.02, G9.62+0.20, T-Tau/Sb, AFGL 2591, NGC 7538, IRAS 06061+2151, and AFGL 2789. Error bars 
show statistical uncertainties obtained with the MCMC analyses. 



Table 4. Correlation coefficients between Galactic parameters for result 4 (49 sources) 





Rq 


UsiC 


»src 


W src 


a 


Rq 


1.000 -0.164 


0.096 


-0.397 


0.065 


0.494 


n 


1.000 


-0.576 


-0.052 


-0.006 


-0.048 






1.000 


0.027 


0.023 


0.002 


^rc 






1.000 


-0.031 


-0.148 


Ware 








1.000 


0.022 


a 










1.000 



Table 5. MCMC results of the Galactic parameters for different model parameters 



result 


#of 


i?o 




a or (a , b Q ) c 




V SIC 


Wsrc 


Vq 


xVd.o.f. 


ID 


source 


(kpc) 


(km/s/kpc) 




(km s" 1 ) 


(km s" 1 ) 


(km s" 1 ) 


(kms^ 1 ) 




10 


52 


8.27±0.42 


29.98±0.64 


+0.05±0.02 


1.1±1.3 


-14.8±1.2 


-0.6±1.0 


5.25* 


494/150 


11 


52 


8.25±0.40 


29.15±0.63 


+0.05±0.02 


1.1±1.4 


-8.1±1.2 


-0.5±1.0 


12.0* 


497/150 


12 


52 


8.25±0.40 


28.10±0.66 


+0.05±0.02 


1.3±1.4 


0.0* 


-0.5±1.0 


19.9±1.2 


499/150 


13 


49 


7.84±0.40 


30.61±0.75 


+0.00±0.02 


0.6±1.4 


-13.1±1.2 


-1.9±1.0 


5.25* 


353/141 


14 


49 


7.82±0.41 


29.60±0.74 


+0.00±0.02 


0.8±1.4 


-6.3±1.2 


-1.9±1.1 


12.0* 


353/141 


15 


49 


7.77±0.38 


28.89±0.71 


+0.00±0.02 


0.7±1.4 


0.0* 


-1.9±1.1 


18.2±1.2 


357/141 


16 


46 


7.94±0.42 


30.59±0.76 


+0.01±0.02 


1.2±1.4 


-13.2±1.3 


-1.3±1.1 


5.25* 


312/132 


17 


46 


7.93±0.40 


29.74±0.74 


+0.01±0.02 


1.3±1.4 


-6.5±1.3 


-1.3±1.1 


12.0* 


314/132 


18 


46 


7.90±0.38 


28.86±0.76 


+0.01±0.02 


1.3±1.4 


0.0* 


-1.3±1.1 


18.3±1.2 


315/132 


19 


49 a 


7.68±0.28 


30.61±0.53 


-0.00±0.02 


0.5±1.0 


-12.8±0.9 


-1.8±0.8 


5.25* 


603/141 


20 


49 h 


7.97±0.51 


30.50±0.97 


+0.01±0.03 


1.0±2.0 


-13.3±1.7 


-1.9±1.5 


5.25* 


192/141 


21 


49 


7.72±0.41 


30.58±0.71 


-0.1±0.7 


0.6±1.4 


-13.2±1.3 


-1.9±1.0 


5.25* 


361/140 










0.1±0.2 












22 


49 


7.70±0.40 


29.71±0.71 


-0.1±0.7 


0.7±1.4 


-6.5±1.3 


-1.9±1.1 


12.0* 


367/140 










0.1±0.2 












23 


49 


7.70±0.40 


28.91±0.71 


-0.1±0.7 


0.6±1.4 


0.0* 


-1.9±1.1 


18.3±1.3 


364/140 










0.1±0.2 













* : fixed parameter 

a : with a v = 5 km s _1 instead of 7 km s _1 . 
b : with cr„ = 10 km s _1 instead of 7 km s -1 . 

c : through models 10 to 20, a (non-dimensional RC index), and for models 21-23, upper raw corresponds to gq = dO/dR 
(km/s/kpc), and lower raw corresponds to b = (d 2 <d/dR 2 )/2 (km/s/kpc 2 ), respectively (values defined at R = R ) 
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some correlation coefficients exceed 0.2, indicating relatively 
large correlation between parameters. Part of such correlation 
can be seen, for instance, in figures 2a, 2c and 2f, where cases 
with higher Rq for results 1-3 provide higher a and lower V STC , 
but note that the results presented here are determined includ- 
ing the effect of such correlation. Of course, the strongest cor- 
relation exists between Vq and V^ rc (which is not shown in 
table 4), which prevents us from determining these two param- 
eters at the same time. 

As is described in section 2, the full source list consists of 
52 sources for which results of precise astrometry are avail- 
able. However, in the sample sources, there likely exist "out- 
liers," which have large deviations from the model-predicted 
radial velocity and/or proper motions. To investigate the ef- 
fect of such outliers, we ran MCMC analyses for several sets 
of sample sources by eliminating possible outliers one-by-one. 
The elimination was done based on the conditional probabil- 
ity Pi(Oj|M): a new sample of TV — 1 sources was created by 
eliminating the source with the smallest value of Pj(Oj|M) in 
the MCMC results for the TV-source sample. In table 3 (and 
5) we also provide x 2 an d degrees of freedom (d.o.f.) for each 
model (rather than the likelihood), because they are fairly use- 
ful to see the effect of outliers. 

As we will see later in this section (e.g, table 5), the pa- 
rameters related to the Galactic rotation velocity such as Qo, 
V SIC and Vq are tightly correlated with each other, and thus 
will be treated in further detail in next subsection (by solving 
for V SIC or Vq). Other parameters such as Rq, a, U STC , and 
W SIC are independent of the choice of model parameters, and 
thus the results listed in table 3 give strong constraints on these 
parameters. In figures 2, we show resultant model parameters 
with the number of sample sources N STC varying from 52 to 44: 
Ro, Qq, a, U SIC , V SIC , Wstc, and \ 2 versus the number of the 
sources N SIC . For reference, the sources with small conditional 
probability are also labeled in table 1. The effects of "out- 
liers" can been seen in figures 2. For instance, in figure 2d the 
value x 2 rapidly increases beyond 50, and x 2 of the full sample 
(52 objects) is nearly 500 while it is around 350 for 49 objects 
(see also results 1 and 4 in table 3). Therefore, it is likely that 
at least these three sources (IRAS 20126+4104, G48. 6 1+0.02, 
G9. 62+0. 20) are indeed outliers, and these three sources may 
be better to be excluded from the sample to obtain the best es- 
timates of the Galactic parameters. In fact, previous studies re- 
ported large peculiar motions for these sources: G48.61+0.02 
was found to be associated with supernova W51C (Nagayama 
et al. 2011), and thus the peculiar motion is likely to be orig- 
inated from a supernova explosion. Sanna et al.(2009) sug- 
gested that large deviation (~ 40 km s" 1 ) of G9. 62+0. 02 from 
the Galactic rotation may be due to the bar potential in the in- 
ner Galaxy. The cause of relatively large non-circular motion 
of IRAS 20126+4104 is unknown, but Moscadelli et al.(201 1) 
already pointed out an existence of that at ~20 km s _1 level 
when a standard rotation curve is considered. 

Due to the effect of these outliers, the best-estimate values of 
some parameters such as Rq, a, etc. indeed vary with decreas- 
ing TVgrc from 52 to 49 (as in figures 2). On the other hand, re- 
ducing N stc further (below 48) do not alter the best-estimates 
of those parameters as much. Therefore, one can expect that 
the sample with ~48 sources do not contain "outliers." 
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Although at least the top-three sources in terms of low con- 
ditional probability are highly likely to be outliers, here we 
do not simply eliminate these outliers from our samples, but 
we take a rather conservative approach by taking both sam- 
ples with/without outliers into account and consider the effects 
of outliers as possible systematic error in our analyses. For 
instance, in figure 2a, the Galaxy center distance Rq is de- 
termined with an statistical uncertainty (T stat of 0.4 kpc (e.g., 
for N src ~ 48). On the other hand, the best estimated of Rq 
for 7V src = 50 and 52 differs from that for TV src = 48 slightly 
beyond the statistical error bar. This clearly demonstrates the 
systematic effect of outliers, and we regard this scatter as the 
systematic error in our estimates of the Galactic parameters. 
In practice, we obtain the best estimate of the parameters by 
taking the mean of the results with the sample source number 
varying from 44 to 52, and also calculate its standard deviation 
as the systematic error er sys . Then we evaluate the final un- 
certainty by summing up in quadrature both typical statistical 
error (taken as the mean of statistical error bars for N SIC rang- 
ing from 44 to 52) and systematic error cr sys . Note that here 
the selection of the smallest number of the sample sources to 
be used is rather arbitrary, but in the present analyses we adopt 
44 for the minimum number of the sources. In fact lowering 
the number of the sample further would not alter the resultant 
parameters drastically, but the uncertainty would be slightly in- 
creasing due to statistical effect. 

The best estimates of Galactic parameters obtained in that 
way are also plotted in each panel of figures 2 by dotted lines 
showing the possible range. For the best estimate of Rq, we 
obtained 

Ro = 8.05 ±0.45 kpc. (8) 

This is 5%-level determination of the most fundamental 
Galactic constant Rq. As we will see in discussion later, this 
is consistent with other recent studies of Rq such as using mo- 
tions of stars around the Galaxy center or variable stars, but 
we note that our result is independent of others and provides 
another constraint on the distance to the Galaxy center. 
The rotation curve index a is also tightly constrained to be 

a = 0.022 ±0.029. (9) 

The index is within the error bar, confirming that the rotation 
curve of the Galaxy is basically flat in the regions where our 
samples are distributed (Galacto-centric radius of 4 to 13 kpc). 
The use of the full sample gives a marginally positive value of 
a of 0.05, but we note that a = 0.05 results in a rise of rotation 
velocity only by ±3 % between Galacto-centric radius between 
5 and 15 kpc. Hence, even in the maximum case of a, Galactic 
rotation curve flatness is ±3% level in those regions. 

The U and W components of the source non-circular mo- 
tions are also well constrained here, and we obtained 

U SIC = 0.95 ±1.46 km s" 1 , (10) 

Were = -1.41 ±1.19 km s" 1 . (11) 

The mean inward motion U SIC is obtained to be within the 
error bar. Also, the mean perpendicular motion W STC is found 
to be practically within the error bars, being independent of 
the models. This is rather natural when one considers the sym- 
metry of Galactic disk along the Galactic latitude. 
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These parameters do not depend on the various models con- 
sidered in the next section, and we take them as our final re- 
sults for these parameters. Regarding other parameters, from 
figures 2 we also obtained Hq = 30.43 ± 0.80 km s _1 kpc -1 , 
V STC = — 13.60 ± 1.49 km s -1 , but as will be seen in the next 
subsection, we have to consider other systematic effects to ob- 
tain the final results. 

5.2. Results with changing model parameters 

In this subsection we present the MCMC results of the 
Galactic parameters by choosing different model parameter 
sets and also different rotation curves. Table 5 summarizes the 
results, again presenting results of the best estimates of the pa- 
rameters, as those in table 3. In table 5 results 10-12 are a 
family of the best estimates for the full sample of 52 sources, 
with different treatment of V^ rc and Vq: in result 10 we adopt 
^=5.25 km _1 (Dehnen & Binney 1998, note that result 10 is 
identical to result 1 in table 3), while in result 1 1 Vq=12.0 km 
~ 1 (Schonrich et al. 2010) is adopted. In both results 10 and 11, 
V^ rc is considered as a parameter to be solved for. As seen in 
results 10 and 11, the SFRs lag velocities (negative V^rc) differ 
by nearly the same amount with the change in adopted Vq, in- 
dicating strong correlation between V SIC and Vq. On the other 
hand, in result 12, instead we fixed Kj r c=0 km s _1 (no SFRs 
lag) and solved for the Solar motion Vq. In other words, in re- 
sult 12, the apparent slow rotation of SFRs is totally attributed 
to the Solar peculiar motion. Results 13-15, and 16-18 are 
other families of the results similar to results 10-12, but for 
different numbers of the sample sources: results 13-15 are for 
the sample of 49 sources, and results 16-18 are for the sample 
of 46 sources (note that results 4 and 13, and results 7 and 16 
are identical to each other, respectively). The elimination of 
sources from the full sample is conducted based on the condi- 
tional probability P;(0;|M) , as is already described in 5.1. 
These results with different number of sources are again pre- 
sented to evaluate the effect of samples with/without outliers. 

Results 19 and 20 are to investigate the dependence of the re- 
sults on the adopted value of 1 -dimensional velocity dispersion 
of SFRs (a v ): throughout the paper, a v =l km s _1 is adopted, 
but only for results 19 and 20, this is altered to 5 or 10 km s _1 , 
respectively (Vq=5.25 km s" 1 is adopted for both results 19 
and 20). Results 19 and 20 should be compared with results 
13, because they are based on the same sample and parameters 
with only difference in a v . As can be seen in table 5, the results 
19, 20 and 13 are fairly consistent with each other. The only 
difference is the size of statistical error bars and the value of 
\ 2 , as is expected in case of using different a v : result 19 gives 
X 2 of 603 against d.o.f. of 141, and result 20 gives x 2 of 192. 
They are to be compared with result 13 provides x 2 of 353. In 
terms of reduced x 2 (= X 2 /d-0.f), a v of 10 km s _1 seems too 
large, and our choice of a v = 7 km s _1 seems more reasonable. 
The reduced x 2 may indicate that a v = 5 km s _1 is reasonbale, 
but here we conservatively adopt a v = 7 km s _1 because the 
use of too small a v could lead to underestimation of statistical 
error. Note that the adopted value of a v = 7 km s _1 is similar 
to those in the previous studies: for instance, Reid et al. (2009b) 
used 7 km s _1 for a v , and McMillan & Binney (2010) solved 
a v by setting it as a parameter to be estimated, and obtained a v 
from 7 to 10 km s _1 depending on their models. 




V® (km/s) 

Fig. 3. Correlation between Vq and V sr c ■ Open circles are those for 
results 10-12 in table 5, rilled circles with error bars for results 13-15, 
and open triangles for results 16-18, respectively. Thin line corre- 
sponds to the best fit relation of V ar c = Vq — 19(±2) km s — 1 , with 
dotted lines showing the error bar range. 

The last three results (21-23) in table 5 are those for us- 
ing different rotation curve models. Reid et al. (2009b) and 
McMillan & Binney (2010) suggested that the Galactic param- 
eters also depend on choice of rotation curve shape. In most 
cases in table 3 and 5 we used a power-law rotation curve. In 
order to see how results would be altered with a different rota- 
tion curve, we adopt the 2nd-order-polynomial rotation curve 
in equation (2) for results 21-23. This rotation curve has two 
parameters, while a power-law rotation curve has only one pa- 
rameter. However, as can be seen in table 5, the use of a differ- 
ent rotation curve does not change our results drastically. This 
may indicate that the shape rotation curve determined with the 
current sample is basically well-represented by a simple func- 
tion such as the one-parameter power-low model. 

From table 5, one can see that the parameters Rq, a, U SIC 
and W SYC are not strongly dependent on the choice of Vq or 
V^ rc . For instance, from results 11 to 20, Ro ranges from 7.7 
to 8.3 kpc, but basically this is the sample effect: in fact, the 
best estimate of Rq remains unchanged within a family of re- 
sults based on the same sample (such as results 10-12, 13- 
15 etc.). The variation of Ro for results 10-23 thus originates 
from the difference in the number of sources (and hence inclu- 
sion/exclusion of outliers). The best estimate of Rq obtained 
in 5.1, which is i?o=8.05±0.45 kpc including systematic error, 
fully covers the range of Ro in table 5. Therefore, as we de- 
scribed in the previous subsection, we can consider the value 
of Rq obtained in the previous subsection as the final value for 
the current study. In addition to Rq, the same applies to a, U sxc 
and W STC , and the best estimates obtained in section 5.1 cover 
the full ranges of these parameters in table 5. 

On the other hand, in table 5, there are strong correlations 
between the Solar peculiar motion Vq, the mean source mo- 



8 



Honma et al. 



[Vol. , 



33 



32 



u 31 



30 



a 29 



28 



27 



o 



10 15 

V© (km/s) 



o 



20 



u 

a 

on 

£ 

G 



25 



Fig. 4. Correlation between and S1q . Open circles are those for 
results 10-12 in table 5, rilled circles with error bars for results 13-15, 
and open triangles for results 16-18, respectively. 

tion in the rotational direction V src , and the angular velocity of 
the LSR Hq. One can see this correlation by comparing the 
results for the same sample but with different Vq or V^c (e.g., 
results 10-12). To clearly show this situation, in figure 3 we 
show the relation between Vq and V^ rc for the results in table 
5. The open circles in figure 3 are for results 10-12, filled cir- 
cles (with error bars) for results 13-15, and open triangles for 
results 16-18, respectively. As seen in figure 3, the relation be- 
tween Vq and V SIC can be approximated by a linear function, 
with V SIC = Vq — 19 km s _1 (thin line in figure 3). When one 
compares \ 2 values with the three results based on the same 
sample (e.g., models 10-12, 13-15, ...), there is little difference 
in x 2 , indicating that our results do not show any preference for 
a particular choice of V^ rc and Vq . Therefore, at this stage, the 
only thing we can constrain is that these two quantities must 
satisfy the following relation 

-l 



y src = y -19(±2)kms- 



(12) 



The allowed region of (Vq, V aIC ) with the above error bar 
(±2 km s _1 ) is also shown in figure 3 with dotted lines. If 
Vq = 5.25 km s _1 is adopted to this relation, the star-forming 
regions lag behind the Galactic rotation by ~ 14 km s _1 , con- 
firming the earlier suggestion by Reid et al. (2009b). On the 
other hand, an upward modification of Vq by Schonrich et al. 
(2010) with Vq = 12.0 km s _1 reduces the lag down to ~ 7 km 
s _1 level. Only Vq of ~ 20 km s" 1 can lead to no lag case, 
but this solar motion is larger than the classical IAU value of 
the Solar motion, with \V\ = 20 km s" 1 toward (18'\ 30°) 
in B1900 (e.g., Kerr & Lynden-Bell 1986), corresponding to 
Vq = 15 km s _1 . Hence, this should be taken to be an extreme 
case, and probably reality would lie among these cases (e.g., 
up-ward modified Vq with some rotation lag). The tight cor- 
relation between Vq and V sr c can be simply attributed to the 
fact that the observed sources motions are relative to that of the 
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Fig. 5. Angular rotation velocity of the Sun, Qq, obtained for results 
10-23 in table 5 (filled circles with error bars). Horizontal axis is the 
result ID, running from 10 to 23. Dotted lines show the possible range 
of Qq for our best estimate, Qq = 31.09 ± 0.78 km s -1 kpc — . 

Table 6. Additional Galactic parameters 



model 


Vq 


Vq - V 8IC 


n Q 


e 


ID 


km s 1 


km s _1 


km s 1 kpc 1 


km s _1 


10 


5.25 


20.1±1.2 


30.61±0.64 


247.9±13.7 


11 


12.0 


20.1±1.2 


30.60±0.63 


240.5±12.8 


12 


* 


19.9±1.2 


30.50±0.69 


231.8±12.5 


13 


5.25 


18.4±1.2 


31.28±0.75 


240.0±13.6 


14 


12.0 


18.3±1.2 


31.13±0.74 


231.5±13.4 


15 


* 


18.2±1.2 


31.23±0.74 


224.5±12.3 


16 


5.25 


18.5±1.3 


31.25±0.76 


242.9±14.2 


17 


12.0 


18.5±1.3 


31.25±0.74 


235.8±13.3 


18 


* 


18.3±1.2 


31.18±0.78 


228.0±12.5 


19 


5.25 


18.0±0.9 


31.29±0.53 


235. 1± 9.5 


20 


5.25 


18.5±1.7 


31.16±0.97 


243.1±17.4 


21 


5.25 


18.5±1.3 


31.26±0.71 


236.1±13.7 


22 


12.0 


18.5±1.3 


31.27±0.71 


228.8±13.1 


23 


* 


18.3±1.3 


31.29±0.74 


222.6±12.8 



* denotes that Vq is solved for by setting V STC = km s . 
The resultant Vq is identical to the value in the right column of 
Vq — V SYC . 0o is calculated using the relevant values of Rq and 
Slo in table 5. 

Solar system and hence the effect of the Solar motion and lag 
of star-forming regions are still degenerated. 

The correlation between V^c and Vq also affects the angular 
rotation velocity of the LSR Sl - For instance, f2 simply de- 
creases in results 10-12 as Vq increases, and the same trends 
are found in results 13-15 and 16-18 in table 5. Figure 4 il- 
lustrates this correlation between Vq and firj, with the same 
symbols as in figure 3. Here again Vq and £Iq show a linear 



No. ] 

Table 7. Summary of Galactic parameters obtained in the present pa- 
per. EiTor bars indicated here include systematic errors. 



Ro 


8.05±0.45 


kpc 






31.09±0.78 


km s" 


1 kpc 


Q 


0.022±0.029 








0.95±1.46 


km s" 


■ l 


v SIC 


Vq - 19 (±2) 


km s" 


■ l 


w src 


-1.41 ±1.19 


km s" 


■ l 



trend. Considering this situation, it is useful to introduce a fol- 
lowing new Galactic parameter 

"o = «o + ^, (13) 
no 

which is the angular velocity of the Sun. Table 6 summarizes 
the derived Q and other parameters such as Vq — V src and 
Go, based on the results in table 5. Also, figure 5 shows f2 
for results 10-23 in table 5. As seen in figure 5, the resultant 
f2 shows only weak dependence on the results (number of the 
source, choice of Vq etc.), and thus is readily constrained well 
in the present study. As was done for Rq determination in the 
previous subsection, from figure 5, f2 is determined as 

n & = 31.09 ±0.78 km s -1 kpc" 1 . (14) 

Here again the error bar includes systematics due to the ef- 
fect of outliers. This is consistent with the direct measure- 
ment of Sgr A* motion along the Galactic plane carried out by 
Reid & Brunthaler (2004), who obtained m = 6.379 ± 0.024 
mas yr" 1 , or the angular rotation velocity fi of 30.24±0.11 
km s" 1 kpc" 1 . Note that while f2 is well constrained, in table 
6 9o still remains relatively uncertain. In fact, in table 6, the 
full possible range of 0rj is 223 - 248 km -1 . Hence the de- 
termination of precise Vq and Oo is still an important issue for 
future observations. In table 7, we summarize final values and 
relations for Galactic parameters in the present study. The er- 
ror bars in table 7 include both statistical and systematic error, 
as already discussed in this section. 

6. Galactic parameter comparisons 

6.1. R 

In the previous section, we have obtained the Galaxy center 
distance to be Ro = 8.05 ± 0.45 kpc, with the error bar in- 
cluding systematics due to possible "outliers". This result is 
fully consistent with recent studies of Ro- For instance, the 
review paper by Reid (1993) summarized the previous stud- 
ies of Ro and provided Ro = 8.0 ± 0.5 kpc. Observations of 
stellar orbits around Sgr A* also provided the distance to the 
Galaxy center as R = 8.4 ± 0.4 kpc (Ghez et al. 2008) and 
Ro = 8.33 ± 0.35 kpc (Gillessen et al. 2009). Earlier analy- 
sis of Galactic structure based on maser astrometry (Reid et al. 
2009b) also reported Rq = 8.4 ± 0.6 kpc. Period-Luminosity 
relation for Mira also suggests R — 8 . 24 ± .08 1 sta t ± .42 1 sys 
kpc (Matsunaga et al. 2009). All these studies are consistent 
within error bars and thus the result in the present paper con- 
firms that Ro is between 8 and 8.5 kpc. 
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6.2. Qq and Qq 

The most direct determination of fi is the apparent proper 
motion of Sgr A* along the Galactic plane: assuming that the 
Sgr A* is at rest at the center of the Galaxy, its proper motion 
directly reflects the sum of the Galactic rotation at the LSR 
and the Solar motion with respect to the LSR. As already men- 
tioned in the previous section, our result of fi = 31.09 ± 0.78 
km s" 1 kpc" 1 is consistent with the Sgr A* motion obtained 
by Reid & Brunthaler (2004), which was ft of 30.24±0.11 
km s" 1 kpc -1 . Note that our astrometry data does not include 
the proper motion of Sgr A* itself, and hence these two stud- 
ies are totally independent. This consistency assures that the 
Sgr A* is indeed at rest at the Galaxy center, with an upper 
limit of its in-plane motion less than ~ 6 km s" 1 . This limit 
is still larger than the upper limit for the proper motion in b- 
direction (perpendicular to the plane), which is found to be ~ 1 
km s" 1 . However, the upper limit of Sgr A*'s in-plane motion 
is significantly smaller than 0o itself (which ranges 223-248 
km s" 1 in the present paper), indicating that Sgr A* is truly at 
the dynamical center of the Galactic rotation. 

Miyamoto & Zhu (1998) analyzed HIPPARCOS data of OB 
stars and Cepheid variables, and derived Q = 32.7 ± 1.4 
km s" 1 kpc" 1 from the OB stars and Q = 30.4 ± 1.4 km s" 1 
kpc" 1 from the Cepheids, respectively. Note that here the re- 
sults by Miyamoto & Zhu (1998) are converted to f2 by us- 
ing their own determination of the Solar motion, which was 
Vq ~ 15 km s . Our result of f2 lies in the middle of 
these two results and consistent with both of them. Earlier 
studies of Galactic parameters based on VLBI astrometry also 
resulted in similar values. For instance, Reid et al. (2009b) 
obtained £7 =30.3±0.9 km s" 1 kpc" 1 based on Vq = 5.25 
km s" 1 , yielding fi = 30.9 ± 0.9 km s" 1 kpc" 1 . Reanalysis 
of the same astrometry data also resulted in similar range of 
O , with 30.9 to 32.5 km s" 1 kpc" 1 (McMillan & Binney 
2010) and 32.9±1.22 km s -1 kpc" 1 (Bobylev & Bajkova 
2010). These results also support that f2 slightly larger 
than 30 km s" 1 kpc" 1 . On the other hand, if one adopts 
IAU recommendations of Galactic constants, one would obtain 
r2 =(22O+15)/8.5=27.6 km s" 1 kpc -1 . Therefore, the present 
paper as well as recent other studies on f2 strongly suggest 
that the Galactic angular velocity should require an upward 
modification by 10% level. Because 0o is tightly linked to f2o, 
the same applies to the rotation velocity of the LSR Oo: while 
the resultant values of Oo are dependent of Vq as well as Rq, 
our results of ©o range between 223 km s" 1 and 248 km s" 1 , 
being higher than the IAU recommendation of 9o=220 km s -1 . 

6.3. Source peculiar motion and Vq 

As we have found in the previous section, the U and W com- 
ponents of the mean source peculiar motion are independent of 
Vq and hence well constrained. As summarized in table 7, the 
mean perpendicular motion W src (= — 1.41 ± 1.19 km s" 1 ) can 
be regarded as at 1.2-a level, which is naturally expected 
from the equilibrium of Galactic disk. The U component is 
also within the error bar, and we can safely conclude that 
the U component is negligible when compared to the Galactic 
rotation velocity, which is an order of ~ 240 km s -1 . 

On the other hand, as we have seen in the previous section, 
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the V component of the source mean motion is tightly corre- 
lated with Vq , and thus exact value cannot be fixed based solely 
on the VLBI astrometry data. The only constraint we can ob- 
tain is the linear relation between Vq and V SIC given in equa- 
tion (12). If we adopt Vq = 5.25 km s _1 of Dehnen & Binney 
(1998), then V^ rc becomes — 14±2 km s _1 . Therefore, our re- 
sult confirms the findings by Reid et al. (2009b) that the star- 
forming regions lag behind the Galactic rotation. It is remark- 
able that this confirmation is done based on much larger num- 
ber of sources compared to that in Reid et al. (2009b). However, 
as we have already seen, models without V SIC also provide a 
reasonable fit to the data: \ 2 only slightly vary with changing 
Vq or solving for it with assuming V 3IC = km s _1 (see, e.g., 
results 12, 15, and 18 in table 5). Therefore, at this stage we 
cannot conclude whether or not a Galaxy model with the rota- 
tion lag of the star-forming regions is better than a model with 
larger Vq and no lag. However, the strong relation in equa- 
tion [12] will provide a good base of future investigations on 
this problem: once one of the parameters (Vq or V STC ) is deter- 
mined based on other observations etc., one can immediately 
find the conclusion on the (non)existence of the SFRs lag. 

6.4. Rotation Curve 

While a few Galactic parameters still have a considerable 
uncertainty depending on the Solar motion etc., the shape of the 
rotation curve itself is well determined based on the analyses in 
the previous section. Figure 6 shows an example of the rotation 
curve, based on result 1 . Here we adopted from table 3 Galactic 
parameters of Ro = 8.27 kpc and Oo = 29.98 km s _1 kpc -1 
(or 0o=248 km s _1 ), respectively. As seen in the figure, the 
rotation curve can be basically regarded as "flat" between 4 
kpc and 15 kpc, being consistent with the rotation curves of 
other galaxies with rotation velocities comparable to that of the 
Milky Way Galaxy (e.g., Sofue & Rubin 2001). 

The flat rotation curve is also consistent with previous stud- 
ies of the Galactic rotation curves (e.g., Clemens et al. 1985). 
On the other hand, Honma & Sofue (1997) showed that the 
shape of HI rotation curve is strongly dependent on ©o and 
for a flat rotation curve, their preferred value of 0o is around 
200 km s _1 . This Oo is inconsistent with the results of the 
VLBI astrometry here, which suggested Oo much higher than 
the IAU standard value of 220 km s ~ 1 . For solving this discrep- 
ancy, again the Solar motion would be a clue, because both of 
the HI rotation curve and Galaxy parameters based on VLBI 
astrometry are significantly affected by the Solar motion V©. 

The flatness of the rotation curve strongly constrained here 
can be used to estimate the distribution of dark matter in the 
Galaxy. In fact, in order to sustain this flat rotation curve, a 
considerable amount of dark matter is required. For instance, 
suppose that the Galaxy's disk is an exponential disk with a 
scale length of h = 3 kpc and maximum disk rotation velocity 
equal to Oo, and the mass-to-light ratio is constant through the 
disk. Then, as is discussed in Honma et al. (2007), the ratio 
of disk mass to the total mass within 13 kpc is around 0.7, 
and hence 30% of the total mass should be in the form of dark 
matter. 

Note that the data points in figure 6 mostly lie below the 
best-fit rotation curve. This is obviously the effect of the lag 
of the star-forming regions because in figure 6 we adopt the 



Solar motion of 5.25 km s _1 that makes the lag largest among 
the results in table 5. If larger Vq is adopted, then the lag will 
be reduced and the data points should scatter around the model 
curve. Other interesting feature in figure 6 is that the devia- 
tion of the sources from the rotation curve could be systematic 
(such as lower rotation velocity at Galacto-centric radii of 5 
kpc and 9 kpc). Although the number of the sample sources is 
not enough to discuss this at the present, this may be related 
to non-axisymmetric structure of the Galaxy (such as a bar and 
spiral arms), and so this will be interesting target for future 
studies with more sources. 
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